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Accurate modeling of cables is important to study the behavior of high 
frequency disturbances in power converter systems. This paper reviews and 
compares two popular methodologies to model energy cables - an improved 
per unit length parameters based model and a Laplace SPICE element based 
model. The two models presented take into account the frequency 
dependence of the parameters of the cable. A ladder network is used for this 
purpose in the per unit length parameters based model. The Laplace SPICE 
element model is generated from a rational function approximation for the 
admittance parameters that are frequency dependent. The rational function 
approximation is obtained using a well known vector fitting algorithm. The 
time and frequency domain solutions of a two wire energy cable, obtained 
from the two models, agree well. 
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1. INTRODUCTION 

Power electronic systems in industrial applications often use long cables. Semiconductor devices 
used in such systems generate pulses with short rise and fall times and thus give rise to high frequency 
components in the currents and voltages [1]. The high frequency common mode and differential mode 
currents in cables produce emissions that can cause unwanted electromagnetic interference (EMI) [2]. 
Therefore, an accurate model of the cables that describes their high frequency behavior is necessary to 
compute the currents flowing in them and thereby predict any emissions produced. The power electronic 
semiconductor devices for drives or other applications are often modelled in SPICE [3]. The high frequency 
behavioral models of motors with the models of the power electronic system and cables are simulated in 
SPICE to predict the electromagnetic interference [4]. Hence, a high frequency circuit model of cable in 
SPICE is necessary. The signals flowing on long cables of power electronic systems are pulse width 
modulated (PWM). Such signals can be used as excitations in SPICE to test the cable model when the cable 
is terminated in any load. 

Many existing models are available on standard SPICE simulators for two wire cables. One model 
based on transmission line equations simulates a lossless uniform transmission line. This model describes the 
cable by its characteristic impedance and signal propagation delay. Since this model does not include the 
cable losses, it is not suitable for use in applications to predict EMI when long cables are used. Another 
simple SPICE model is based on cascaded sections of pi, T or L-shaped equivalent circuits of per unit length 
R, L, G and C parameters distributed along the cable. The RLGC model ignores the variation of the per unit 
length resistance and inductance with frequency due to the skin effect and proximity effect. Similarly, the 
variation of the permittivity and loss tangent of the dielectric, which results in increased loss at higher 
frequencies, is not accounted for in this model. By ignoring those parametric variations, inaccurate time and 
frequency domain solutions are obtained from simulations. Several methods have been proposed to correct 
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this by using improved models for the cable. Such models account for the variation of the parameters with 
frequency. One method has been the use of a ladder network to model the frequency dependence of losses 
[5]. Another method has been to use a rational function approximation of the admittance matrix of the cable 
[6-12]. These two methods have been presented and compared in this paper. 


2. THEORY 

Accurate analysis of EMI in power converter systems for drives requires high frequency circuit 
models of all the components in the system. A commercial tool, ANSYS-Simplorer that integrates multiple 
modelling technologies has been used to simulate the power electronic converter system [13]. Further this 
work used, finite element analysis based ANSYS-Maxwell to obtain the high frequency behavioral model of 
induction motor and Matlab-Simulink to simulate its control and study faults caused by different drive’s 
components. The high frequency behavioral models of all the components in the power electronic systems 
including the cables can also be obtained from different software platforms and simulated in SPICE to predict 
EMI [4]. For this approach, two methods to obtain a SPICE cable model is reviewed and presented below. 

The electrical properties of a transmission line are characterized by the per unit length parameters R, 
L, G and C. The circuit theory-based definitions of the parameters are related to the stored electric and 
magnetic energies, as well as the power loss and obtained from a field theory-based analysis of the 
transmission line [14]. The transmission lines parameters of some common lines such as coaxial, two wire 
and parallel plate are determined at low frequency using analytical techniques. The lumped parameter SPICE 
model of the transmission line is obtained by using a number of L-shaped sections of the RLGC parameters 
distributed along the length of the cable. If n sections (L-shaped) are used per unit length as shown in Figure. 
1, the lumped element values in each section are R/n, L/n. G/n and C/n. The first and last series branches are 
half sections with lumped elements R/(2n) and L/(2n) to keep the line symmetric. There should be enough 
sections in the model to ensure that each section represents the length of the line that is a small fraction of the 
wavelength (A/20) at the highest frequency of the signal flowing on the line. Depending on the rise time of 
the time domain signal, the highest frequency component on the line is estimated as 035/(rise time). Some 
designers use the estimate of l/(rise time) for the highest frequency component. In practice, by choosing the 
number of sections to be large enough in the model, the system behavior with a long cable is predicted in the 
time and frequency domains. However, this model is still insufficient to describe the effects of cable length 
on signal flow, as the per unit length parameters are frequency dependent. In order to determine the 
parameter values at high frequencies, numerical techniques, such as the finite element method, are used [15]. 
Their variations with frequency are later accounted in the SPICE model. One method to account for the 
frequency dependence of the parameters is the use of a ladder network in the model [5]. The losses due to 
skin effect, causing the cable behavior to be frequency dependent, are modelled by the RL ladder network. 
This ladder network is used in place of a single R and L in the series branch of Figure 1. The frequency 
dependent variation of the real and imaginary components of permittivity of the dielectric material of the 
cable is modelled with an RC ladder network. This ladder network is used in place of a single G and C in the 
shunt branch of Figure 1. 



Figure 1. Transmission Line Model with n Sections 
of Per Unit Length Parameters R, L, G and C 
Distributed on Unit Length of the Cable 



Figure 2. Cross-section of the Two Wire Cable: 
Diameter of the Wires D 2 = 1.38mm, Diameter of 
the PVC Covering of Wires Di = 2.88mm, 
Diameter of the Rubber Shield Over the two 
Wires D 3 = 8.7mm, Relative Permittivity of 
Rubber = 2.3, Conductivity of the 
Wires =45.94MS/m 
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Another method to account for frequency dependence of per unit length parameters is by directly 
using frequency dependent elements in the SPICE model. To do this, Laplace SPICE elements are used to 
represent the frequency dependent admittance parameter matrix of the long cable. A vector fitting algorithm 
approximates the admittance parameter matrix of the cable as a function of frequency [9] with a rational 
function in the Laplace domain. A Laplace SPICE element is then used to represent the rational function 
approximation. A model that characterizes the transmission line with the frequency dependent effects is thus 
obtained. The results from the ladder network SPICE model and Laplace element SPICE model are compared 
here. The ladder network model for a two wire cable, as presented in [5], has been calculated here for the 
purpose of comparing the two methods. 

2.1. Ladder Network Model for a Two Wire Cable 

The cross section of a 1 m two wire unshielded energy cable is shown in Figure 2. Using a finite 
element method-based electromagnetic solver tool (ANSYS-HFSS), the RLGC parameters of this cable were 
obtained over the frequency range lOKHz - 100MHz. These parameters can be also determined with a vector 
network analyzer. The lm long two wire cable was modelled with number of sections, n = 32. The time 
domain pulse with a rise time of 10ns was used to test the model (highest frequency component of 100MHz). 
With 32 section per meter of cable, the length of each section was much lower than A/20 at 100MHz. The 
variations of resistance and inductance with frequency (series branch of Figure 1) are represented by an RL 
ladder network shown in Figure 3. A least square curve fitting algorithm was used to determine the elements 
of the ladder network. This algorithm minimized the difference in equivalent impedances of the ladder 
network and the series branch of the cable in the frequency range. The variations of conductance and 
capacitance parameters (shunt branch of Figure 1) with frequency are represented by the RC ladder network 
shown in Figure 4. Using the same least square curve fitting algorithm, the difference in admittances of the 
ladder network and the shunt branch were minimized. A SPICE sub circuits netlist for Figure 3 and Figure 4 
are written with the determined ladder network element values. The SPICE model for the cable is created by 
including 32 sections of the sub circuits in the netlist. 

The model has been validated by comparing the input impedances of the lm-long cable terminated 
with either a short or open circuit. The experimental results for these configurations from [5] have been 
compared with the model obtained in this work as shown in Figure 5. Although the recalculated values of the 
elements of the ladder network vary slightly from the original work, the simulated impedances are in good 
agreement with the measurements. 

The time domain response of the cable to a 40V pulse with a rise and fall time of 10ns is shown in 
Figure 6. The applied input 40V pulse was identical to the pulse in [5] and the transient response simulation 
from this work agrees well with measurements [5]. This solution is further compared with the model obtained 
from the admittance parameter matrix of the cable. 



Figure 3. RL ladder network. For the lm cable of 
Figure. 2 with 32 sections, Rl=1.5e-1Q, 
R2=9.9021e-2Q, R3=7.5125e-2Q, R4=4.25e-2Q, 
R5=6.9246e-3Q, Ll=5.0002e-13H, L2=6.0733e- 
11H, L3=6.0733e-11H, L4=3.4392e-10H, 
L5=1.7158e-08H. 


Figure 4. RC ladder network. For the lm cable of 
Figure. 2 with 32 sections, Rl=6.5199e8Q, 
2=4.561 le6Q, R3=7.0987e6Q, R4=1.3797e7Q, 
R5=4.2357e9Q, Cl=0.54018e-15F, C2=1.3475e- 
15F, C3=2.9342e-15F, C4=1.9026e-15F, 
C5=1.65e-12F. 
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2.2 Laplace SPICE element model for a two wire cable 

The admittance matrix of a network relates the current to the voltages at each of its ports [14]. A two wire 
cable of any length forms a two port network. For a two port network, the equation below shows the relationship 
between current and voltage at the input (port 1) and output (port 2). 

7i Y lt Y 12 1 [Pil 
U 2 J“^2i Y 22 \[v 2 \ 

Since the cable is a reciprocal network, Y 12 = Y 21 . The admittance parameters matrix of the cable 
can be measured using a vector network analyzer over a desired frequency range. This requires placing 
appropriate connectors between the cable and the network analyzer and impedance termination of 50 Q at the 
end of cable. The admittance parameters matrix of the two wire cable can also be obtained from a finite 
element analysis software. For a lm long cable with the cross-section shown in Figure 2, the admittance 
parameter matrix has been obtained from a finite element solver (ANSYS-HFSS) over a frequency range of 
1 OKHz-100MHz. 




Figure 5. Input Impedance of lm Cable Terminated Figure 6. Transient Response of lm Cable 

in Short Circuit (SC) and Open Circuit (OC) Terminated in a 750 Q Load, Connected to 50 Q 

Source and an Input Pulse of 40V (Starts at t = 30 
ns with Rise Time of 10ns, Duration 150ns, 
Fall Time of 10ns) 


The admittance parameters obtained are a function of frequency and therefore can be expressed as a 
function in the Laplace domain (s = i * 2nf ). A rational function approximation of the admittance matrix 
T(s) is determined by using the vector fitting algorithm [9]. The admittance parameters, or in other words, 
the elements of the matrix T(s) of the network are expressed in pole residue form by this algorithm and 
are given as: 


iV 

Y (s) = V — l,mn + D 

trunks) 


Where the N poles s ± to s N are either real or complex conjugate pairs, their corresponding residues 
R t are either real or complex conjugate pairs, and D is real. The vector fitting algorithm expresses all the 
parameters in the admittance matrix with the same set of poles and different residues unless two of the 
parameters in the admittance matrix are equal. 

For a circuit simulation, the rational function approximation should satisfy stability, causality, and 
passivity criteria. The vector fitting algorithm [9] creates an approximation that is stable and causal. The 
algorithm presented in [16] ensures passivity by checking whether the eigenvalues of Re[Y(s)] matrix are 
positive for all frequencies. If the criterion is not satisfied, it uses a perturbation technique to reformulate the 
rational function approximation. 

A two port admittance matrix can be represented by an equivalent network in pi topology [14]. The 
admittance parameters form the branches in the pi topology network. Traditionally, standard network 
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elements represent the admittance parameters. When rational function approximations of admittance 
parameters are in the Laplace domain, the branches of the pi topology can be replaced with these functions as 
shown in Figure 7. Each of the terms in the rational function approximation can be represented as Laplace 
elements in SPICE. By using a pi topology, there are no SPICE nodes inside the two ports of the cable. 
However, there are as many parallel Laplace SPICE elements between the nodes as the number of poles 
(corresponding to the terms in summation) used in approximating the admittance parameter. For example, the 
Laplace SPICE element of Y mn corresponding to, say, the first (i = 1) pole s ± (assumed to be real) is given 
as: 

GXXX m n LAPLACE {V(m,n)} {R 1>mn /(s - s t )} 

For a complex conjugate pole pair, say, s 2 and its conjugate, the Laplace SPICE element is given as: 

GXXX m n LAPLACE {V(m,n)} {(A * 5 + E)/(s 2 + B * 5 + C) J 
where the values of B, C, A, and E are: 

B = -2*Re(s 2 ) 

C = |s 2 | 2 
A = -B 

E = -2* Re(R 2 s 2 ) 

The constant D mn is represented in SPICE as: 

GXXX mn LAPLACE {V(m,n)} {D mn } 

All the poles of the rational function and are expressed as G elements (Voltage dependent 
current source) in the SPICE netlist. Similarly, all branches in the pi topology are expressed as G elements to 
create SPICE model of the cable. For the lm-long two wire cable, a rational function approximation with as 
few as 20 poles is sufficient to model the cable. The computed and the measured [5] values of impedance 
when terminated in short and open circuits, as shown in Figure 2, are found to agree well. The time domain 
response of the cable to a 40V pulse with rise and fall times of 10ns is shown to agree well with 
measurements [5] and the ladder network representation. 

The computation was repeated for a 10m long cable of the same cross section. Both models were 
tested in time domain by applying a 40V pulse of 10ns rise time. The results obtained from the two models, 
as shown in Figure 8, agree well. 



i -Yi 2 (s) 2 



Figure 7. Two Port Network Model of two Wire Figure 8. Transient Response of 10m Cable 

Cable with Admittance Parameters in Laplace Terminated in a 750 Ohm Load, Connected to 50 Ohm 
Domain Source and an Input Pulse of 40V (Starts at t=40ns 

with Rise Time of 10ns, Duration 190ns, Fall Time of 
10ns) 


3. CONCLUSION 

Two SPICE models have been obtained to simulate a system that uses a long cable and switching 
pulses with short rise times. The LTspice simulator was used to simulate the system. From the simulation of a 
long cable (lm and 10m), it is observed that both models present solutions with good accuracy. The ladder 
network model is composed of several internal nodes in the cable, while the Laplace SPICE element model 
requires no internal nodes to represent the cable. The number of internal nodes increases with the length of 
the cable in the ladder network model. The SPICE circuit model netlist for the cable using Laplace elements 
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has fewer lines. However, the computation times to perform a time domain simulation using the Laplace 
SPICE element model are longer than those of the ladder network model. This is due to the fact that while 
performing transient analysis in SPICE, the impulse response of the Laplace elements are calculated and 
convolved with input signal. If the model arrived at using Laplace elements is not passive, numerical 
inaccuracy results. If steps are taken to avoid this problem when synthesizing the model, accurate transient 
analysis can, however, be realized. Although either of the methods used to obtain a SPICE model for the 
cable is appropriate in terms of accuracy, the use of vector fitting algorithm to obtain a model is found to be 
elegant and likely to be more suitable for multiconductor cables. 
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